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Optimal forwarding ratio on dynamical networks with heteroge- 
neous mobility 
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Abstract -As the discovery of non-Poissonian statistics of human mobility trajectories, more 
attention has been paid to understanding the role of these patterns in different dynamics. In this 
study, we first introduce the heterogeneous mobility of mobile agents into dynamical networks, and 
then investigate the forwarding strategy on the heterogeneous dynamical networks. We find that 
the faster speed and the higher proportion of high-speed agents can enhance the network through- 
put and reduce the mean traveling time in the case of random forwarding. A hierarchical structure 
in the dependence of high-speed is observed: the network throughput remains unchanged in small 
and large high-speed value. It is interesting to find that the slightly preferential forwarding to 
high-speed agents can maximize the network capacity. Through theoretical analysis and numerical 
simulations, we show that the optimal forwarding ratio stems from local structural heterogeneity 
of low-speed agents. 



Introduction. — The notion of mobile agent is really 
important in the study of many dynamical systems, such 
as collective motion in biology [1], information transmis- 
sion on wireless ad- hoc network [3] , and traffic congestion 
in city transportation system [3] . In all of these dynamical 
systems, each mobile agent is free to move independently 
in any direction, and locally interacts with current nearby 
agents. Take wireless ad- hoc network for example [2]. The 
continuous movement of mobile devices would frequently 
change their links to other devices due to limited forward- 
ing range. Thus, these dynamical systems can be treated 
primitively as dynamical networks, in which connections 
between nodes always change temporarily Network the- 
ory provides therefore a natural framework to study the 
emergence of structural and dynamical properties of these 
systems [Ml]. 

Physicists recently became increasingly fascinated by 
the dynamics on dynamical networks, such as epidemic 
spreading of mobile individuals [9l413j , synchronization of 
mobile oscillators [T4"Ml6) . and consensus of communicat- 
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ing agents [T7| . Especially for packet forwarding of mobile 
devices, how to enhance transmission efficiency has be- 
come a crucial problem due to the rapid development of 
wireless communication technology. Considering costs and 
technology limitations, more efforts should be made to im- 
prove the routing strategy [2fT8H23] . Yang et al. first stud- 
ied transportation dynamics on dynamical networks [24], 
and compared random routing with greedy routing }25| . 

However, the previous researches have only focused on 
the case of identical mobile agents, which can not reflect 
the diversity of mobile individuals. Since the discovery 
of non-Poissonian statistics of human behaviors such as 
human interaction activities |26j and mobility trajecto- 
ries |27H28j . more and more scientists have been paying 
attention to the role of these patterns in different dynam- 
ics. Most recent research results showed that both time 
and space activities have significant impacts on spreading 
dynamics [291 - 136] . In this letter, we introduce the het- 
erogeneous mobility of mobile agents into dynamical net- 
works, and then propose an optimal forwarding strategy 
on the dynamical networks with heterogeneous mobility. 
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Model. — We consider N agents performing random 
walk in an area of L x L with periodic boundary conditions. 
Initially, the agents are randomly put on the area. Let 
Vi(t) be the speed and 9i(t) be the angle with respect to 
the cc-axis characterizing the direction of the velocity of 
agent i, and Xi(t) and yi(t) be the coordinates of the agent 
i and time t. The coordinates and the velocity of agent i 
at time t + 1 are updated according to 



Xi(t+1) = Xi(t) + Vi(t)cos9i(t), 

Vi(t + 1) = yi (t)+Vi(t)ainSi(t), 

Vi(t+1) = Vi, 

0i(t + l) = fc(t), 



(1) 



where £j (t) is a random variable chosen uniformly between 
the interval [— 7r, n] after the coordinates arc updated. To 
account for heterogeneous mobilities of the mobile agents 
(27] [28], we assume that the agents only take on either 
a low speed vi or a high speed v^. The fraction of high 
speed agents is /, i.e., the system is characterized by a 
distribution of speeds D{v) among the agents represented 

by 

D(v) = (l-f)S(v-v l ) + fS(v-v h ), (2) 

where S(x) is the Dirac <5-function. The speed of an agent, 
once assigned initially, does not vary with time. There- 
fore, an agent moves with a constant speed but in random 
directions. 

To consider the packet routing efficiency in a network of 
mobile agents, let R packets be generated at a time step. 
Each of these packets is generated randomly at some agent 
i and placed at the end of the queue of agent i, with a des- 
tination that is also randomly chosen among the agents in 
the system. To forward the packet to its destination, each 
agent is capable of forwarding at most C packets in its 
queue to selected neighbors every time step. The delivery 
follows a First-In-First-Out policy. As the agents are mo- 
bile and in view of the connectivity of electronic devices, 
we regard the agents that arc instantaneously within an 
area of radius r of an agent i to be the neighbors, i.e., the 
neighbors of agent i are those with dij < r, where dij is the 
distance between agent i and agent j. The radius r thus 
sets the range over which a packet could be forwarded in 
a time step. In this point of view, a mobile-agent network 
corresponds to a dynamical network, with the links con- 
stantly broken and established as the agents move. The 
algorithm for forwarding a packet goes as follows. If the 
destination of a packet is found among the neighbors, the 
packet will be delivered to its destination and removed 
immediately. If not, an agent decides whether to forward 
it to a high-speed neighbor with a probability p or to a 
low-speed neighbor with a probability 1 — p. After mak- 
ing a decision, a neighbor of the chosen type is randomly 
picked and then the forwarded packet is put at the end of 
the queue of the selected neighbor. If there is no neighbor 
of the chosen type, the packet will not be delivered. The 
process is repeated for each packet to be forwarded and 



every agent carries out the algorithm in every time step. 
As there is a fraction / of high-speed agents in the system, 
the case of p = f roughly equal to the case of randomly 
choosing a neighbor to forward a packet regardless of its 
type, and the cases of p < f (p > /) correspond roughly to 
forwarding a packet to low-speed (high-speed) neighbors 
preferentially. 

Let S(t) be the total number of packets in the system 
at time t. In the congested phase, there will be more 
and more packets accumulated in the system; while in 
the balanced phase, the number of delivered and removed 
packets balances that of packets generated. These two 
phases can be characterized by [23] 



— lim 

R t-too 



(S(t + At) - S(t)) 
At 



(3) 



which plays the role of the order parameter in that r\ > 
in the congested phase and r\ = in the balanced phase. 
For a given value of C, the system is in the balanced 
(congested) phase for packet generation rates R < R c 
(R > R c ). Asa figure of merit for packet delivery, a 
higher R c corresponds to a better algorithm. 

Simulation Results. — First we investigate the high- 
speed agents how to influence network throughput R c in 
the case of random forwarding, i.e., p = f. Fig. 1(a) 
depicts the phase transition (from free flow state to con- 
gestion state) for different Vh values when / = 0.4, and R c 
is found to depend onw/,. A hierarchical structure in the 
dependence is shown in Fig. 1(b): R c remains unchanged 
in small or large u/, value, which is consistent with the 
results in Ref. [21] ■ Besides, Figs. 2 (a) and (b) show 
that R c increases with /, which denotes more high-speed 
agents can enhance the network throughput. Then the 
question is: Why does the speed Vh and the proportion 
/ of high-speed agents influence the network throughput? 
We will discuss it later (see Eq. (fl~6|0 . 

Next, we think of the case of preferential forwarding, 
i.e., p^/. As shown in Fig. 3(a), it is very interesting to 
note that there is a maximum R c value when p is slightly 
greater than /, e.g., R m sw 15 at p m as 0.48 when / = 0.4 
and Rm « 22 at p m w 0.88 when / = 0.8. As p m > f 
corresponds to the preferential forwarding to high-speed 
agents in Fig. 3(b), this optimal phenomenon means that 
the slightly preferential forwarding can enhance network 
throughput to some extent. For example, R m » 22(p s=s 
0.88) > R c ks 16(p = 0.8) when / = 0.8. Moreover, 
Fig. [4] shows that faster speed and the higher proportion 
of high-speed agents can reduce the mean traveling time 
(T) (defined as the average number of hops for all data 
packets from their sources to destinations) in the case of 
random forwarding. 

Theory. — In order to understand the above optimal 
phenomenon, here we present a mean field analysis of net- 
work throughput. As all agents are randomly distributed 
on the planar space at any time t due to the random walk 
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Fig. 1: (Color online) The network throughput as a func- 
tion of high speed value for the case of random forwarding, 
(a) Order parameter rj(R) versus R for different Vh values 
when / = 0.4. r)(R) is obtained by averaging over 10 3 
time steps after disregarding 3.5 x 10 3 initial steps as tran- 
sients, (b) The network throughput R c versus Vh when 
/ = 0.2, 0.4, 0.8. The lines are the theoretical predictions 
from Eq. (fT0| where (T), kf max and k^ nax are obtained 
by numerical simulations. The parameters are chosen as 
N = 10 3 ,L = 10, r = I, vi = 0.001, C = l,p = f. The 
results are obtained by averaging over 10 independent re- 
alizations. 



Fig. 2: (Color online) The network throughput as a func- 
tion of the proportion of high-speed agents for the case of 
random forwarding, (a) Order parameter rj(R) versus R 
for different / values when vh — 2. i](R) is obtained by 
averaging over 10 3 time steps after disregarding 3.5 x 10 3 
initial steps as transients, (b) The network throughput 
R c versus / when Vh = 0.01,0.1,1,5. The lines are the 
theoretical predictions from Eq. (|TU|) where (T), kf max 
and fc£ max are obtained by numerical simulations. The 
parameters are chosen as TV = 10 ,L = 10, r = l,vi = 
0.001, C = l,p = f. The results are obtained by averag- 
ing over 10 independent realizations. 



of agents, the degree distribution will be approximated 
by [37] 



P(k) 



~W(fc) fc 
k\ 



(4) 



where (k) = Nitr 2 /L 2 . 

In the balance phase, according to the forwarding strat- 
egy, the evolution of the packet numbers of low-speed 
agent I with ki(t) and high-speed agent h with kh(t) can 
be written as 



dni(t) 
dt 

dn h {t) 
dt 



(1 - p)A ljnj (t) 
(!-/)%(*) ' 



N 



3=1 



(5) 



where the sum runs over all nodes of the network, and Aij 
is the element of the adjacency matrix, in which Aij = 1 
when dij < r, otherwise Aij = 0. At large time t, all 
agents will reach the balance of traffic flow, and we have 



N 



(l-p)A ljnj (t) 
(!-/)%(*) ' 



rihit) 



N 

E 



pA hj nj(t) 



(6) 



Considering Nf high-speed agents randomly distributed 
on the planar space, we know that an agent with k(t) 
has fk(t) high-speed neighbors and (1 — f)k(t) low-speed 
neighbors. We suppose ni(t) = B\ki{t) and nh(t) = 
B2kh{t), and then Eq. ((6]) is transformed into 



nt(t) 
n h (t) 



{ ^ft[(l-f)B 1 + f B 2 ], 
pkh(t) . 



f 



■[(l-f)B 1 + fB 2 \, 



(7) 



where the occupation packet number is independent of 
degree correlations |37| . From the above, we know that 
Bi/B 2 = [(1 - p)/]/[p(l - /)], and Eq. © is also written 



ni{t) = 
n h {t) = 



B{l- P )k,{t) 

1-/ ' 

Bpk h (t) 

f ' 



(8) 



where B is a constant. In the balance phase, there are 
S = J2i n i(t) + ^2h n h{t) — R(T) packets on the network 
at each time step [55]. Submitting Eqs. (U) and ([5J into 
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Fig. 3: (Color online) The network throughput as a func- 
tion of the forwarding rate when / is fixed, (a) The net- 
work throughput R c versus p when / = 0.2, 0.4, 0.8. The 
lines are the theoretical predictions from Eq. (|16|) where 
(T)(C = oo ) is obtained by numerical simulations (see 
Fig. [3]) and kf max is given by Eq. (fT4"|). A difference be- 
tween the theoretical predictions and the numerical re- 
sults is originated from the real (T) rea i > (T)(C = oo) at 
the critical point, (b) The optimal forwarding p m versus 
/. The red solid line is the theoretical prediction from 
Eq. (|T7|) where kf max is given by Eq. (fT4")l . and the blue 
dotted line corresponds to p = f. The parameters are cho- 
sen as N = 10 3 ,L = 10, r = l,v t = 0.001,% = 2,C = 1. 
The results are obtained by averaging over 10 independent 
realizations. 



Fig. 4: (Color online) The mean traveling time (T) as a 
function of the forwarding rate p in the free flow state 
with C = oo when / = 0.2, 0.4, 0.8. The inset shows (T)p 
as a function of the forwarding rate p. The parameters 
are chosen as N = 10 3 ,i = 10, r = l,vt = 0.001,% = 2. 
The results are obtained by averaging over 10 independent 
realizations. 



space in 0((T)) steps due to their fast moving. According 
to the law of large number, we have 



lim 



1 



o«D)->oo 0((T)) 



(11) 



this self-consistent relationship, and we obtain 



B 



R(T) 
N(k)' 



(9) 



At the critical point, the dynamical balance of traf- 
fic flow happens in a relatively long time with the same 
order as 0((T)). The node possessing the maximum 
occupation packet number must meet the relationship 
Et° i{T))n raa X (t) < CO«T» in 0((T)) steps. According 
to Eqs. © and ©, we have 



R ( T )P k h, m ax 

N(k)f 
R{T)(l-p)kl max 
N(k)(l-f) 



<a 



(10) 



where k\ max and kf max are the maximum effective de- 
gree of high-speed agents and low-speed agents respec- 
tively, and the effective degree of node i is defined as 
= o((T)) St ki(t). 

Owing to the heterogeneous mobility, kf max of low- 
speed agents is different from that of high-speed agents. 
All high-speed agents can span a broad area of planar 



The maximum effective degree of high-speed agents fc£ max 
is approximativcly equal to (k). On the other hand, low- 
speed agents are taken as motionless nodes in this stage 
because of its very small speed. For this reason, the effec- 
tive degree of a low-speed agent is divided into two parts: 
f(k) high-speed neighbors and k[ settled low-speed neigh- 
bors. The number of settled low-speed neighbors follows 
a Poisson distribution 



z-^){k[f 

k[\ 



(12) 



where (ki) = (1 - ])N ixr 2 j 1? . Thus, the distribution of 



kf follows 



P(kf) = 



-( k 'i){ki) (kf ' f{k)) 



(13) 



(*?-/(*»! ' 
The natural cut-off effective degree of low-speed agents 
can be calculated by [33] 



N(l - f) / P(kf)dkf 
'fcf_„ 



1. 



(14) 



Submitting Eq. (TTTj) into Eq. (ITU1) , the network through- 
put is given by 



R c = Min{R h c: R l c }, 



(15) 
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and 



Rl- 



NCf 
OV 

N(k)C(l - f) 
{T){l-P)kf 



(16) 



where Min{} represents the less of the two. If p = f, 
Eq. (fTB")) will go back to the case of random forwarding. 
When Vh is much smaller than r, the dynamic network is 
approximatively taken as a fixed structure in a compara- 
tively long time, so the network throughput changes very 
little. As Vh increase to the value with the same order as r, 
the shorter mean traveling time, which is due to the longer 
forwarding distance in each hop, results in the increase of 
network throughput. For large values of Vh, the network 
throughput remains unchanged because of the fixed mean 
traveling time and maximum degree |24j . For this reason, 
a hierarchical structure in the dependence of high-speed is 
observed in Fig.[TJb). In this case, (T) (due to more high- 
speed agents) and kf max (from Eq. (|T4"))) decrease with /. 
Thus, R c = R l c increases with /, which is validated by the 
results in Fig. [21b). From Eq. (fTU)) , the theoretical pre- 
dictions are basically consistent with the numerical results 
in Figs. [Tf b) and[2jb). When Vh is not very large (e.g., 
Vh =1), there is a difference between the theoretical pre- 
dictions and the numerical results due to the invalidation 
of effective degree hypothesis. 

When / is fixed, (T) decreases with p because the pref- 
erential forwarding to high-speed agents can increase the 
distance which packets are forwarded in one time step (see 
numerical results in Fig. [4]) . From Eq. (fTB]). we know that 
decreases with p because of the increase of (T)p (con- 
firmed by numerical results in the inset of Fig. [4]), while 
R l c increases with p because of the decrease of (T)(l — p). 
There exits a maximum network throughput R m when 
R l c = R^ (see Fig. [3] (a)), and the corresponding optimal 
p m is thus obtained by 



Pn 



fK 



(i-.f)(k) + fkf 



(17) 



Interestingly, p m only depends on the structure of the dy- 
namical network. Owing to the Poisson characteristic of 



effective degree of low-speed agents, i.e., kf 



> {k), Pn 



will be slightly greater than /, which means that the opti- 
mal forwarding rate is the slightly preferential forwarding 
to high-speed agents. As shown in Fig. 0(b), the theoreti- 
cal predictions are in good agreement with the simulation 
results. 

Conclusion and Discussion. In conclusions, we 
have investigated the forwarding strategy on the dynami- 
cal networks with heterogeneous mobility. First, we have 
shown the faster speed and the higher proportion of high- 
speed agents can enhance the network throughput and 
reduce the mean traveling time in the case of random for- 
warding. A hierarchical structure in the dependence of 
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Fig. 5: (Color online) Effective degree distribution for two 
kinds of agents when / = 0.2 (a) and / — 0.8 (b). The 
parameters are chosen as N = 10 3 ,L = 10,?' = l,vi = 
0.001,% = 2,0«T)) = 50. The results arc obtained by 
averaging over 10 2 independent realizations. The lines are 
the theoretical predictions from Eq. (TlT 



high-speed is observed: the network throughput remains 
unchanged in small and large high-speed value. Second, it 
would be interesting to know that the slightly preferential 
forwarding to high-speed agents can maximize the network 
capacity. By combining theoretical analysis with numeri- 
cal simulations, we have obtained the optimal forwarding 
rate, and found this optimal phenomenon stems from local 
structural heterogeneity of low-speed agents. This work 
provides us further understanding and new perspective in 
the effect of human dynamics on the forwarding strategy 
on dynamical networks, and thus may help to optimize the 
delivering strategy on heterogeneous dynamical networks. 



Ming Tang would like to thank Pakming Hui for stimu- 
lating discussions. This work is supported by the NNSF of 
China (Grants No. 11105025), China Postdoctoral Science 
Foundation (Grant No. 20110491705), the Specialized Re- 
search Fund for the Doctoral Program of Higher Education 
(Grant No. 20110185120021), China Postdoctoral Science 
Special Foundation (Grant No. 2012T50711), and the 
Fundamental Research Funds for the Central Universities 
(Grant No. ZYGX2011J056). 

REFERENCES 

[1] VlCSEK T. and Zafeiris A., Phys. Rep., 517 (2012) 71. 
[2] Abolhasan M., Wysocki T. and Dutkiewicz E., Ad 

Hoc Networks, 2 (2004) 1. 
[3] Barthelemy M., Phys. Rep., 499 (2011) 1. 
[4] Albert R. and Barabasi A.-L., Rev. Mod. Phys., 74 

(2002) 47. 



p-5 



Yu Gan et al. 



[5] Boccaletti S., Latora V., Moreno Y., Chavez M. 
and Hwang D.-U., Phys. Rep., 424 (2006) 175. 

[6] Dorogovtsev S. N., Goltsev A. V. and Mendes J. F. 
F., Rev. Mod. Phys., 80 (2008) 1275. 

[7] Barrat A., Barthelemy M. and Vespignani A., Dy- 
namical Processes on Complex Networks (Cambridge Uni- 
versity Press, Cambridge, UK) 2008. 

[8] Newman M. E. J., Networks: An Introduction (Oxford 
University Press, Oxford, UK) 2010. 

[9] Frasca M., Buscarino A., Rizzo A., Fortuna L. and 
Boccaletti S., Phys. Rev. E, 74 (2006) 036110. 

[10] Buscarino A., Fortuna L., Frasca M. and Latora 
V., Europhys. Lett, 82 (2008) 38002. 

[11] Peruani F. and Sibona G. J., Phys. Rev. Lett., 100 

(2008) 168103. 

[12] Li X., CAO L. and Cao G. F., Eur. Phys. J. B, 75 (2010) 
319. 

[13] Yang H.-X., Wang W.-X., Lai Y.-C. and Bing-Hong 

Wang, Europhys. Lett, 98 (2012) 68003. 
[14] Frasca M., Buscarino A., Rizzo A., Fortuna L. and 

Boccaletti S., Phys. Rev. Lett, 100 (2008) 044102. 
[15] Fujiwara N., Kurths J. and Diaz-Guilera A., Phys. 

Rev. E, 83 (2011) 025101. 
[16] Frasca M., Buscarino A., Rizzo A. and Fortuna L., 

Phys. Rev. Lett, 108 (2012) 204102. 
[17] Baronchelli A. and Di'az-Guilera A., Phys. Rev. E, 

85 (2012) 016113. 
[18] Echenique P., Gomez-Gardenes J. and Moreno Y., 

Phys. Rev. E, 70 (2004) 056105. 
[19] Echenique P., Gomez-Gardenes J. and Moreno Y., 

Europhys. Lett, 71 (2005) 325. 
[20] Yan G., Zhou T., Hu B., Fu Z. -Q. and Wang B. -H., 

Phys. Rev. E, 73 (2006) 046108. 
[21] Zhang H., Liu Z., Tang M. and Hui P. M., Phys. Lett. 

A, 364 (2007) 177. 
[22] Tang M., Liu Z., Liang X. and Hui P. M., Phys. Rev. 

E, 80 (2009) 026114. 
[23] Tang M. and Zhou T., Phys. Rev. E, 84 (2011) 026116. 
[24] Yang H.-X., Wang W.-X., Xie Y.-B., Lai Y.-C. and 

Wang B.-H., Phys. Rev. E, 83 (2011) 016102. 
[25] Yang H.-X., Wang W.-X., Lai Y.-C. and Wang B.-H., 

larXiv:1112.4773V l (2011). 
[26] Barabasi A.-L., Nature, 435 (2005) 207. 
[27] Brockmann D., Hufnagel L. and Geisel T., Nature, 

439 (2006) 462. 
[28] Gonzalez M. C, Hidalgo C. A. and Barabasi A.-L., 

Nature, 453 (2008) 779. 
[29] Vazquez A., Racz B., Lukacs A. and Barabasi A.-L., 

Phys. Rev. Lett, 98 (2007) 158702. 
[30] Iribarren J. L. and MORO E., Phys. Rev. Lett., 103 

(2009) 038702. 

[31] Tang M., Liu L. and Liu Z., Phys. Rev. E, 79 (2009) 
016108. 

[32] TANG M., Liu Z. and Li B., Europhys. Lett, 87 (2009) 
18005. 

[33] Wang P., Gonzalez M., Hidalgo C. and Barabasi 

A.-L., Science, 324 (2009) 1071. 
[34] Balcan D. and Vespignani A., Nat. Phys., 7 (2011) 581. 
[35] Belik V., Geisel T. and Brockmann D., Phys. Rev. X, 

1 (2011) 011001. 
[36] Zhao Z.-D., Liu Y. and Tang M., CHAOS, 22 (2012) 

023150. 



[37] Antonioni A. and Tomassini M., larXiv:1207.2573H 
(2012). 

[38] Zhao L., Lai Y.-C, Park K., and Ye N., Phys. Rev. E, 

71 (2005) 026125. 
[39] Boguna M., Pastor-Satorras R. and Vespignani A., 

Eur. Phys. J. B, 38 (2004) 205. 



p-6 



